Induced Pseudoscalar Coupling Constant 



Il-Tong Cheon * and Myung Ki Cheoun ' 
Department of Physics, Yonsei University, Seoul, 120-749, Korea 

(9 November, 1998) 

Abstract 

The recent TRIUMF experiment for yTp — ► nva") gave a surprising result 
that the induced pseudoscalar coupling constant gp was larger than the value 
obtained from fj,~p — > nv^ experiment as much as 44 %. Reexamining contri- 
bution of the axial vector current, we found an additional term to the matrix 
element of Beder and Fearing which was used to extract the gp value from 
the measured photon energy spectrum. This additional term plays a key role 
to restore the reliability of gp(— 0.88m^) = 6.77^(0). 
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In semileptonic weak interaction, the strong force can generally induce four couplings ad- 
ditional to the usual vector and axial vector couplings, i.e. weak magnetic Gm, pseudoscalar 
Gp, scalar Gs and tensor Gt- 

The matrix element of vector and axial vector currents are given as 

(N(p')\V£(0)\N(p)) = u(p')[Gv(q 2 )r + ¥V + <WK^]f n{p) (1) 
(N(p')\A^0)\N(p)) = u(p')[G A (q 2 )r + 2f£V + G T (g 2 KV]7 5 f u(p) , 

where G A (0) = g A (Q), G M (0) = g M (Q), G v (0) = g v (0) and G P (q 2 ) = (^)g P (q 2 ) with the 
nucleon and muon masses, m and m M . r a is the isospin operator. Gs and Gt belong to the 
second class current which has a different G-parity from the first class current, and they are 
assumed to be absent from the muon capture to be discussed in this paper. On the basis of 
the PCAC (Partially Conserved Axial Current), the induced pseudoscalar coupling constant 
is calculated as 

Off) TV) 

g P (-0.88ml) = g - g A (0) = 6.77^(0). (2) 

This value is confirmed by an experiment of the ordinary muon capture (OMC) on a proton, 

However, such kind of determination of gp value induces 25 % uncertainty at least, 
because the momentum transfer is far from the pion pole. It is extremely important to 
obtain a precise value for gp of the weak hadronic current, because it plays a key role in the 
fundamental weak interaction processes. The only way to approach to the pion pole is the 
radiative muon capture (RMC) on a proton, /i~p — > nu^. 

Recently, the TRIUMF group measured the RMC photon energy spectrum and extracted 
a surprising result [0 

g P = g P (-0.88ml)/g A (0) = 9.8 ± 0.7 ± 0.3 . (3) 

It exceeds the value obtained from the OMC as much as 44%. This discrepancy is serious 
because the theoretical value of gp is predicted in a fundamental manner based on the 
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PCAC imposed on the axial vector current and agrees with the OMC value. As long as 
the PCAC is assumed to be creditable, a doubt may be cast on the result of TRIUMF 
experiment. However, the measured photon energy spectrum seems to be reliable in view 
point of their enough experimental experiences in TRIUMF. In order to solve this puzzle, 
one has to reexamine carefully the Beder- Fearing formula |5],|J used to extract the gp value 
from the measured spectrum. 

The chiral perturbation calculation (ChPT) was recently carried out for OMC || and 
well reproduced the PCAC prediction, i.e. gp = 6.77. It is also consistent with the result of 
the heavy baryon chiral perturbation theory |J. Therefore, we believe that the calculation 
based on PCAC is reliable. Since the RMC amplitude is generated merely by the minimal 
coupling procedures from the OMC amplitude 0, the calculation method based on PCAC 
might preserve the confidence even for the RMC. 

The axial current coupled to the electromagnetic field can be derived in an elegant 
manner, i.e. a standard gauge transformation. We start from the ordinary linear cx-model 
given by the Lagrangian 

C = n^d, + g(a + + l -[{d^) 2 + («V) 2 ] - + a 2 ) - ^(vr 2 + <r 2 f . (4) 

Although this Lagrangian is invariant under the SU(2) isovector infinitesimal chiral and 
gauge transformations, v|/ — >• v|/' = (1+27577- ^f — ► iff — ^(l + rys/f- 1), a — > a' = a + fj-n 
and 7? — > 7r = 7r — ffa, it describes only a massless fermion. In order to create the pion 
mass, the chiral symmetry breaking term (a should be included into Cq. Since the vacuum 
expectation value of the a field does not vanish, i.e. (0|er|0) = /„., the a field is shifted as 
a — > a = a — f n and, then, (0|a|0) = is fulfilled. Then, the pion and sigma meson's masses 
can be given as = /i 2 + A/ 2 , m 2 = /i 2 + 3 A/ 2 and ( = /^m 2 , and the nucleon mass is 
m = -gf v . 

The resulting Lagrangian without the symmetry breaking term (a becomes 
C = #[z 7 ^ -m + g(a + ir- 7r 75 )]$ + ^[(^vf) 2 + (d^) 2 } 
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_m^ 2 _ m|„ 2 _ A[A(7f2 + „ 2) + m2] „ _ A (7f2 + _. 2)2 + ^ _ (5) 



By the relation 



exp(^7 5 r • 0) = cos(y-) + i7 5 r • 0sm(y-) = ^-[(T + ?7 5 r • vf] , (6) 

J-7T Jir Jtt Jtt 

where = 0/0, and replacement of a — > <r' = y + f w cos(j-), the Lagrangian can be 
rewritten for g = —m/f n as 

<, = *[<7 l '^ + ffA«=P(T r-0 75 )]^ + ^(^0) 2 , (7) 

where £q = Co — const. 

It is easy to see that the Lagrangian C' holds a global chiral symmetry under the in- 
finitesimal phase transformations ^ — > ^' = (1 + 275 1 • ^)^, 0^0' = — /^^. Even when 
the derivative is replaced by a covariant derivative, i.e. <9^ — > -D^ = <9 M — ?ee M where is 
the photon polarization vector, £q remains in the global chiral symmetry, since higher order 
terms than e 2 and if are ignored and thus, erf yields negligible small quantities. 

For deriving the axial current from the Lagrangian, let us first obtain the extended Euler 
equation by defining the action, 

S = J d\C{U*),D%{x)) , (8) 

where = <9 m ±k m . Under an infinitesimal variation 6<j>i(x), the stationary condition yields 

Since 8{D^(f> i ) = <5(<9 M 0j ± K^fa) = <9^(<50j) ± K^Sfa, the integrand becomes 

dC . , dC . _ . , . , dC . , 
dfc <9(D±0i) *d{D±<l>i) 

^wr {a ^^^ ] ^ +a ' l wkj^- (10) 

Thus, the extended Euler equation is found as 
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Next, we derive the extended Gell-Mann-Levy equations. The variation 5<pi = —r} a (x)F? 
changes the Lagrangian C as C + 5C, where 5C{<pi, D^fa). Then, we have 

d{Dl '(pi) d{D), '(pi) 

where we used the extended Euler equation in (11). Defining ( F?) = , we obtain 

the extended Gell-Mann-Levy equations, 

a ~ d{D^ a ) ' 0*. • (13) 

If = 0, i.e. r/ a is independent of x, eq.(12) reduces to 5C = —(d^A^r/a. Since 5C = 

if the lagrangian holds a global symmetry under the infinitesimal variation, we have d^A% 

= 0. This statement is synonymous with the Noether's theorem. 

Under the local transformations, ^ — > ^' = ^(l + ry 5 | • f}(x)), * — > = (1 + 275 1 • 

— » — ^ — » ^ 

r/(x))\l/, — > = — f^ffix), and — > £> M = <9 M — iee M , the Lagrangian, £ reads as 

— * 1 2 

C = n^D, + 9fnexp(j l5 T • 0)]* - * 7 "75^ ' (^77) + ^OV^ " " /* ■ 

(14) 

For this case, eq.(12) yields 

5C = -f 7"75^ • (^) " ' A(^) • (15) 

Thus, by eq.(13), the axial current and its divergence are given as 

K^-fl^ + U^a , £>< + M£ = 0. (16) 

Here, notice that = + iee^ because = = — iee^. If we add the chiral 

symmetry breaking term —^j-(f 2 to C' , the second equation in (16) yields 

D^A* = -mlM a , (17) 
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which is the extended form of the PCAC. The same equation was also given by Adler J7j. 
As is seen in eq.(16), the axial coupling constant appears to be unity. However, it is well 
known as qa = 1.25. To cure this defect, the following chiral invariant Lagrangian is added 
to C' 0, 



£1 = ^^75^-/^, (18) 
where C\ is determined so as to give qa = 1.25. Then, the axial current becomes 

A£ = g A ^lJ^ + UD^ a , (19) 

where g A = 1 + Cxf%. 

Operating a covariant derivative on eq.(19) and equating it to eq.(17), we find 

- U(dMD» + m 2 )0 a = #[ + »^75^ , (20) 

i.e., in ignoring e 2 -order term, 

U(Ql - 0<t>a = g A [2miM> 75 y* + ief e M 7 M 75y *] ■ (21) 

When the solution of this equation <p a with a definition, 5u/(<? 2 — m^) = -gp/2mm M , is 
substituted into eq.(19), we obtain 

Trip Tftfj, £ 

+!^[*(*>^7»<*>] ' (22) 

where e 2 term is ignored. The fourth term is missing || in the previous calculations Hf§]. 
As will be shown below, it is actually a very important term corresponding to the seagull 
term in the pion photoproduction M. 



It may be possible to introduce another chiral invariant Lagrangian of the form [10 



£ 2 = C 2 ^-(^x^) , (23) 
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but this one does not contribute to the axial current. This fact can be easily seen as shown 
below. When eq.(23) is added to the Lagrangian, the meson field is obtained by means used 
above in the following form instead of eq.(20), 

$ = [1 + flC\B\\g A B + f n C 2 g A (B x B) + f 2 n C 2 2 g A B(B ■ B)] , (24) 

where B = — 2rrm f P gA / -^Jt (^T^Tsf 1 ^)- Since the second term in eq.(24) vanishes, it reduces 
to = g A B. This is exactly the same as that given in eq.(21). 

Following the Fearing's formulation and notation || for diagrams given in Fig.l, one can 
evaluate the relativistic amplitude of RMC on a proton as 

M fi = M a + M b + M c + M d + M e + AM e (25) 



with 



M a = -e a u n T 5 (Q)u p ■ uMl - 7s) ~ I^^ TX > ( 26 ) 



-id 



M 6 = e a L s u n T d (K ^ f, + mp ( 7 a - iih^-kp)^ , 



a a(j A + # + m. 



M c = e a L s u n (-iK n - — fc^) — — r (iT)w p , 

2m„ 2k ■ n 



2Q a + k a g P (K 



2\ 



M d = -e a L 5 u n (— - K ~f 5 )u 

M e = e a L 5 u n {- — a H >y 5 g )u p , 

2m m„ 



9p{K 



2\ 



AM e = -e a Lsu n (^P^Q 6 ^l a )u } 



2mm M 



where 



rid/ \ <5 i SB , s i 9p{q 2 ) ,i 

r (g) = Pv7 + 9/3 + #A7 75 + 9 75 , (27) 

2m m M 

£<5 = ^75(1 — jsju^K = n — p + k and Q = n — p with momenta of neutron, proton 
and photon, n,p and fc, respectively. And m ~ m p ~ m ra . Other constants are taken as 
= 1.0, pa = — 1.25, gM = 3.71, k p = 1.79 and n n = —1.91 0. M e term is originated from 
the third term in eq.(22) and AM e term comes from the fourth term. But the latter, AM e , 



is missing in the paper by Fearing ||f|]. Accordingly, this term was not included in the 
previous procedure of extracting gp value from the experimental photon energy spectrum 
0. The transition rate is given by 

dT RMC aG^m N r i dy kE? 1 ^ ^ 



dk (2n) 2 J-i *W -k(l-y)A s ^ ns 

where a is the fine structure constant, G is the standard weak coupling constant, y = 
k-i>, k max = (W£-m 2 n )/2W , E v = W (k max - k) /[W - k(l- u)}, W = m p + m n - (muon 
binding energy) and |0 M | 2 is the absolute square of muon wave function averaged over the 
proton which is taken SIS db point Coulomb. 

In order to compare to the experimental results, we take the following steps. For liquid 
hydrogen target, muon capture is dominated through the ortho and para p\ip molecular 
states P,P|. Since these molecular states can be attributed to the combinations of hyperfine 



states of up atomic states []TTJ] i.e. single and triplet states, we decompose the statistical 
spin mixture \ J2 sp ins \Mfi\ 2 into such hyperfine states by reducing 4x4 matrix elements to 
2x2 spin matrix elements. At this step, we confirmed that when the AM e term was not 
included, eq.(28) reproduced the curves given in ref. ||. Finally, by exploiting the mixture 
of muonic states relevant in experiments , we calculate the photon energy spectrum. The 
count number of the photons is now expressed as 

N = Z d ^L . (29) 
dk 

Here Z is determined by adjusting the value of dT RMC /dk without AM e term for g P = 
gp(— 0.88m 2 )/g A (0) = 9.8 so as to agree with the best fit curve in ref. 0. With this value 
of Z, we have to examine the case of AM e included. 

Our results are shown in Fig. 2. The solid curve is for the spectrum obtained in ref. [0, 
i.e. the result without AM e term for gp = 9.8. On the other hand, the dotted curve is 
calculated without AM e term for gp = 6.77. This curve is obviously much lower than the 
measured spectrum. When AM e term is taken into account for gp = 6.77, we obtain the 
dashed curve which is very close to the solid curve. The minor discrepancy may be due to 
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the neglect of higher order contribution and other degree of freedom such as A. Our result 
shows that AM e term restores the credit of gp = 6.77. 

The number of RMC photons observed for k > 60MeV is 279 ± 26 and the number of 
those from the solid curve is 299, while our result obtained by integrating the dotted curve 
spectrum is 273. Since the contribution of A degree of freedom is known to be a few percent 
0, it is not included in the present calculation. Vector mesons such as p and u make also 
very small contributions. Higher order terms are pointed out to be insignificant |§. 

The pion field actually interacts in virtual state with the nucleon and therefore the 7rNN 
form factor may be taken into account as an off-shell effect. However, the standard irNN 
form factor fnNNiq 2 ) — (A 2 — m 2 )/(A 2 — q 2 ) with A 2 = m 2 + m 2 participates through an 
effective gp(q 2 ), i.e. gp(q 2 ) = gp(<f)f-KNN{. ( l 2 ) but gives only 4 ~ 5% contribution, because 
the process occurs at low momentum transfer, q 2 = —0.88m 2 . 



Recently, Kirchbach and Riska |12[] proposed a pseudovector form for the pion-induced 
component of the axial current, i.e. the second term is replaced by 9p ^ ^ q^ (7-^)75 instead of 
including the last term in eq.(22). This procedure turns out to give the same result as that 
presented here for the spectrum of yTp —>■ nv^. On the other hand, Jenkins and Manohar 
|T3|] discussed on an effective Lagrangian which contains higher derivative operators and the 



chiral symmetry breaking quark mass matrix. However, in the present theory, this quark 
mass may be absorbed into the pion effective mass. 



As another attempt, the ChPT calculations of RMC have also been carried out |14| , |i5 
but the gp = 6.77 value could not be extracted. Their results are, more or less, the same 
as those of Fearing's calculation ||,|4j], in which the important "Seagull" term was missing. 
Thereby, these calculations may have to be reexamined, taking higher order terms into 
account. It should be noted that the ChPT can satisfy the gauge invariance but it becomes 
obscure if the A degree of freedom is taken into account. 

In the present framework, our calculation shows that gp = 6.77 is reasonable for both 
OMC and RMC on a proton. 
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FIGURE CAPTIONS 



Fig. 1. Standard diagrams describing radiative muon capture on a proton. 

Fig. 2. Photon energy spectrum for triplet states. The solid curve, which is to reproduce 
the experimental data reasonably, were taken from ref. || , i.e. the result without AM e term 
for g P = 9.8. The dotted curve is obtained without AM e term for g P = 6.77. The dashed 
curve is with AM e for gp = 6.77. The dot-dashed curve is calculated with AM e term alone 
for gp = 6.77. (Figure of direct comparison with the experimental data is not presented 
here because of some problems in PS file transform. Please contact to the authors for more 
informations on our results) 
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